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We theoretically propose and experimentally demonstrate a nonclassicality test of single-mode 
field in phase space, which has an analogy with the nonlocality test proposed by Banaszek and 
Wodkiewicz [Phys. Rev. Lett. 82, 2009 (1999)]. Our approach to deriving the classical bound 
draws on the fact that the Wigner function of a coherent state is a product of two independent 
distributions as if the orthogonal quadratures (position and momentum) in phase space behave as 
local realistic variables. Our method detects every pure nonclassical Gaussian state, which can also 
be extended to mixed states. Furthermore, it sets a bound for all Gaussian states and their mixtures, 
thereby providing a criterion to detect a genuine quantum non-Gaussian state. Remarkably, our 
phase-space approach with invariance under Gaussian unitary operations leads to an optimized test 
for a given non-Gaussian state. We experimentally show how this enhanced method can manifest 
quantum non-Gaussianity of a state by simply choosing phase-space points appropriately, which is 
essentially equivalent to implementing a squeezing operation on a given state. 
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Introduction —Nonclassicality of a quantum state is a 
topic of crucial importance that has attracted a lot of 
theoretical and experimental efforts for long. It provides 
not only a profound conceptual framework to distinguish 
quantum phenomena from classical ones, but also an 
essential practical basis for numerous applications, e.g. 
in quantum information processing. An important ap¬ 
proach to studying quantum mechanics in comparison 
with classical mechanics is to adopt a phase-space de¬ 
scription of a quantum state [l| . A wide variety of quan¬ 
tum systems of continuous variables (CVs) can be ad¬ 
dressed in phase space, including quadrature amplitudes 
of light fields, collective spin states of atomic ensembles, 
and motional states of trapped ions, Bose-Einstein con¬ 
densate, or mechanical oscillators, etc. [H. Investigating 
quantum dynamics in phase space has yielded a great 
deal of intuition to quantum-to-classical transition |3| . It 
also offers a powerful tool to treat problems in quantum 
optics 0 and CV quantum informatics Q. 

A clear signature of nonclassicality is the negativity of 
phase-space distribution, which does not exist in classi¬ 
cal probability distributions. However, its demonstration 
typically requires a full reconstruction of Wigner func¬ 
tion and it is of fundamental and practical signifi¬ 
cance to have a sinmler set of measurements manifesting 
nonclassicality [ 3 , 0 , desirably even when the Wigner 
function is non-negative. For instance, every Gaussian 
state possesses a positive-definite Wigner function, which 
restricts a possible set of nonclassicality tests. To man¬ 
ifest the Bell nonlocality, e.g., by employing homodyne 
detections, a Gaussian state must be transformed to a 
non-Gaussian state having a non-positive Wi gne r func¬ 
tion to rule out hidden-variable models I, El. Ba¬ 


naszek and Wodkiewicz (BW) introduced a different sem¬ 
inal approach to manifesting CV nonlocality even with 
a positive-definite Wigner function 11|. Their method 
looks into the particle nature of the field, i.e. the photon- 
number parity directly related to Wigner function, in 
contrast to the wave nature revealed by homodyne detec¬ 
tions. The BW formalism has significantly contributed 
to a deeper understanding of nonclassical correlation, 
which has also been extended to multipartite systems, 
e.g. E- Despite its intellectual impact, the BW 


test was not yet experimentally achieved, furthermore, 
a similar phase-space formalism to demonstrate nonclas¬ 
sicality at a more elementary single-mode level was not 
developed so far. 

Here we propose a nonclassicality test of single-mode 
fields in phase space. Our nonclassicality refers to those 
states that cannot be represented as a mixture of coher¬ 
ent states. We first note that the Wigner function of 
a coherent state is a product of two Gaussian distribu¬ 
tions, leading to a notion of classicality that two orthog¬ 
onal quadratures in phase space behave as local realis¬ 
tic variables. This makes a nonclassicality test, in anal¬ 
ogy with the BW formalism, taking four points in phase 
space. Our test can detect every pure Gaussian nonclas¬ 
sical state. Furthermore, it sets a bound for all Gaussian 
states and their mixtures, thus offering a test of quan¬ 
tum non-Gaussianity—a stronger form of nonclassicality. 
Numerous applications require non-Gaussian resources as 
essential, e.g. universal GV quantum computation [14| . 
entanglement distillation [isjl, quantum error correction 
0 , and GV nonlocality test [l^, etc.. It is thus of cru¬ 
cial importance to distinguish a genuinely quantum non- 
Gaussian state from a mixture of only Gaussian states 
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FIG. 1: (a) Four phase-space points chosen at the vertices of 
rectangle for a nonclassicality test (blue solid) and of paral¬ 
lelogram for a non-Gaussianity test (red solid) (b) An opti¬ 
mal choice of rectangle for testing a squeezed state occurs at 
9 — (p = j {4>: squeezing axis) . 



[13 [3. Note that the latter can also be a non-Gaussian 
state when a finite sum of Gaussian states is taken, giving 
a nonzero value of some known non-Gaussianity measures 
3 . However, it cannot be regarded as a genuine non- 
Gaussian resource. Remarkably, we present a method 
to optimize the non-Gaussianity test using an invariance 
property under Gaussian unitary operations. It enables 
us to detect a broad set of non-Gaussian states, which 
is experimentally demonstrated to illustrate the power of 
our approach. 

Nonclassicality —The Wigner function Wp{q,p) of a 
state p is directly related to the number parity of its dis¬ 
placed state Pa = D{a)pD^{a), where D{a) = “ 

displaces the given state by amplitude a = q + ip Q. 
That is, 


^Wp{q,p) = tr[(-l)“^“pa], (1) 

thus taking a value in [-1,1]. Notably, a coherent state 
I Of) has a symmetric Wigner function that is a product 
of two Gaussian distributions for any pair of orthogo¬ 
nal quadratures, ^W\a){a\{x,y) = 

where {x,y)'^ = TZ{6){q,p)'^ is the coordinate system ro¬ 
tated by an angle 9 from {q,p) in phase space (Fig. 1 
(a)), with ax = Re[ae“*®] and ay = lui[ae~’'^]. The 
Wigner function of a coherent state can thus be re¬ 
garded as a product of two independent random vari¬ 
ables, f Wp(g,p) = ab {0 < a,b < 1). 

Let us define 

1 1 

J[p]^|^^(-l)^'=Wp(x„y,), (2) 

j=0 k—0 


where the four points (xj,yk) form a rectangle oriented 
at angle 9 in phase space (Fig. |T] (a)). We then have 
J[p\ = aobo+aibo+aobi—aibi for a coherent state, which 
is the very structure of the celebrated CHSH inequality 
[3 - Due to the condition 0 < a,6 < 1, we obtain — 1 < 
J <2 2l|. Extending to a mixture of coherent states, 


FIG. 2: (a) Optimal J for a pure Gaussian state ug against 
squeezing ^ = 2tanh2r up to Bo = « 2.32 (gray dashed 

line), (b) Detection of nonclassicality for Gaussian states 
against purity p and squeezing parameter G Our test adopt¬ 
ing four (three) phase-space points detects nonclassicality 
with purity p > « 0.86 {p > 0.5) represented by a purple 

(red) shaded region, (c) Maximal for/|0)(0|-I-(1 —/)|2)(2| 
and (d) e (sinh 7 ^| 0 )( 0 | -I- cosh‘y^\ijj-y){'ijj-y\), using a rect¬ 
angle (blue solid line) and a parallelogram (red dashed line) 
test by optimizing phase-space points for each state. J > 
8/3®^® « 2.32 and J > 2 represent quantum non-Gaussianity 
and nonclassicality, respectively. The Wigner function be¬ 
comes positive for | < / < 1 and 0 < 7 in (c) and (d), 
respectively. 


p = f p(A)|A)(A|, we derive a classicality condition as 


- 1 < J[Pci] < 2 . 


( 3 ) 


We note that ^Wf p(^x)\x){x\ix,y) = Jp{X)a{x\X)b{y\X), 
where the coherent amplitude A plays a role as a hidden 
variable. Violation of inequality (3) indicates a nonclas¬ 
sicality of single-mode fields. 

Our nonclassicality test bears a close connection to the 
BW formalism—a CHSH-type nonlocality test in phase 
space directly addressing two systems [ll|— particularly 
in obtaining the classical upper bound 2. On the other 
hand, our test has a different lower bound -1, of which 
violation also indicates nonclassicality. This necessarily 
requires the negativity of the Wigner function and we 
focus on the violation of the upper bound in this Letter 
[^ . Furthermore, unlike the GHSH inequality allowing 
a maximum 2\/2 (Girelson’s bound) [l^l, the maximum 
value of N in our test approaches 4, which can be at¬ 
tributed to the noncommutativity of the local operators 
q,p and n [22|. 

Gaussian states —We first show how our test detects 
a Gaussian nonclassical state. A single-mode Gaussian 
state can generally be represented as a displaced squeezed 
thermal state. 


(TG = Dia)S{r,(l))ath{n)S\r,(j))D''^{a), 


( 4 ) 
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where S{r, 0) = exp[——is the squeez¬ 
ing operator (r: squeezing strength, (/>: squeezing axis), 
and tTth(h) = (n+i)rt+i l^)(^l is a thermal state (h: 

mean photon number). For a given Gaussian state with 
{a,r, (/), h}, an optimal violation occurs at four vertices 
of a rectangle oriented 45° from the squeezed axis, i.e. 
6 — (j) = j [Fig. 1 (b)] [i^. Intuitively, this choice of 
rectangle makes a sense because our test of nonclassical- 
ity relies on the non-factorizability of Wigner function 
(See the paragraphs below Eq. (1)). The optimal points 
are 


'j o-x — { Vj cty — ( l)-^"*" Zji- 


( 5 ) 


(j = 0,1), with K = ^ 2cosh2r ^nd Zj given in p|. 

In Fig.|2](a), we plot an optimal J[(Jg\ for a pure Gaus¬ 
sian state against the squeezing strength ^ = 2tanh2r, 
whose values are all above 2 demonstrating a success¬ 
ful detection of nonclssicality. To fully identify the de¬ 
tectable regime including mixed Gaussian states, we draw 


a contour plot of J[(tg\ against purity /i = 2 fi+i 
squeezing ^ in Fig.|2](b). Our test detects a mixed Gaus- 

o9/8 

sian state above the purity /i > ~ 0.86 (purple 

shaded region). Remarkably, the detection can be fur¬ 
ther enhanced to ^ > 0.5 (red shaded region) by a sim¬ 
pler test of adopting three phase-space points that has 
no analogue in BW inequality (2^ . 

Among all Gaussian states, we analytically obtain 
the maximal value Bg = ~ 2.32 in the limit 

^ 2 (infinite squeezing), at points (xq,?/ o,a:i,yi) = 

(Kt, —Kt, Sut) with t = [13 • This value co¬ 

incides with the maximal Gaussian bound for two-mode 
[Hj and three-mode [13 nonlocality tests in phase space. 

Genuine quantum non-Gaussianity —A non-Gaussian 
state can take a larger value than the Gaussian bound, 
J > Bg, which becomes a sufficient condition to detect 
quantum non-Gaussianity. As the maximum achieved by 
a Gaussian state ac is Bq, the same bound also applies to 
a mixture of Gaussian states, pmg = Impor¬ 

tantly, under an arbitrary Gaussian unitary operation 
Uq, the Gaussian mixture pmg remains as a Gaussian 
mixture, thus JIIJgPmgUq] < Bq, which enables us to 
construct an optimized test of non-Gaussianity. 
Criterion: Given a state p, if there exists a Uq making 


J[UgpUU >Bg = 8/39/8, 


( 6 ) 


the state p is genuinely quantum non-Gaussian. 

Ug refers to displacement, phase-shift, and squeezing 
operations, however, our test does not require the imple¬ 
mentation of those operations. First, displacement and 
phase-shift do not affect an optimal J[p], thus excluded 
in our consideration. They only shift or rotate the pro¬ 
file of Wigner function while keeping its entire shape, 
thus the optimal points are simply moved along the dis¬ 
tribution but yielding the same maximum 22|. On the 


other hand, under a squeezing operation, using Eq. (1) 
and the identity Syr,(j))D{a)S{r,(j)) = D{a), where 
a = a cosh r -|- sinhr, we obtain W§pgf{x,y) = 

Wp{x,y), with {x,yy^ = S{r, cj)){x,y)'^. The linear 
transformation S{r,(j)) = coshr/ -|- sinh r cos 2(/)cr;; -I- 
sinh r sin 2(/)cra; {ai: Pauli matrix) generally shifts four 
points at the vertices of rectangle to those of parallel¬ 
ogram (Fig. 1 (a)). Therefore, we can essentially imple¬ 
ment the “rectangle” test for the squeezed state SpS"^ 
by simply considering four vertices of the correspond¬ 
ing parallelogram for the given state p, i.e. = 

ij:lk=oi-y^^Wp{x,,y,). 


We illustrate the power of our approach by consider¬ 
ing some examples, compared with other known criteria 
based on the Wi gne r- and the Q-function under energy 
constraint in 2^ [26| , respectively. For a mixture of vac¬ 
uum and two photon states /10) (01 -f (I — /) 12) (21, our cri¬ 
terion detects quantum non-Gaussianity for 0 < / ;< 0.69 
as shown in Fig. 2 (c), whereas the methods of 251 
and [13 detect none and for 0 < / < 0.093, respec¬ 
tively. The Wigner function of a Fock state \n) is given 
by IA|„)<„|(g,p) = f(-l)-e- 29 ^- 2 p^L„( 4 g 2 +V), where 


Ln is the Laguerre polynomial of order n. For the mix¬ 
ture /|0)(0| -I- (1 — /)|2)(2|, the Wigner function be¬ 
comes positive definite for ^ < / < I. Importantly, 
Fig. 2 (c) shows how our detection can be significantly 
enhanced by implementing the parallelogram test (red 
dashed line) in contrast to the rectangle test (blue solid 
line). Similarly, our test detects a mixture of even cat 
state \ijj-f) = A/dq) -f | — 7 )) with coherent amplitude 


7 and vacuum, i.e. e (sinh 7 ^| 0 )( 0 | + cosh 7 ^|^/)-,,)(' 0 .y|) 
for all non-zero 7 (Fig. 2 (d)), whereas the witnesses in 
2 II 23 detect none of these states. 


Experiment —We have performed an experiment for 
the state /|0)(0| -I- (1 — /)|2)(2| using a single mo¬ 
tional mode (X-direction) of a trapped ion in 

3-dimensional harmonic potential with trap frequencies 
(wxiWyjWz) = 27r (2.8,3.2,0.6) MHz. We employ the 
coupling of phonon mode to two electric levels of 
ion in the S' 1/2 ground state manifold, |F = l^mp = 
0 ) = I t) and |F = 0 ,tof = 0 ) = | I) with sepa¬ 
ration whf = (27r) 12.642821 GHz. The control and 
manipulation of phonon mode a is made by the anti- 
Jaynes-Gumming interaction i?ajc = -fh.c. with 

CT+ = I t)(i |i realized by two counter-propagating laser 
beams with beat frequency near whf + wx (22|- Here 
T] = Aky^h/2Mu!x is the Lamb-Dicke parameter, H the 
Rabi frequency of qubit transition, Ak the net wave- 
vector of the Raman laser beams and M the ion mass. 


To test our criterion, we need to (i) prepare an ini¬ 
tial state, (ii) displace it, and (hi) measure the parity, 
sequentially, (i) We prepare the motional ground state 
|n = 0) by Doppler cooling and the standard Raman 
sideband cooling. For the preparation of \n = 2) state, 
we successively apply a 7r-pulse of iJaJG transferring the 
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FIG. 3: (Color online) (a) J for a Fock state |2) by using 
a parallelogram test [Fig. 1 (a)] that essentially implements 
squeezing (r: strength) on a given state. Theoretical predic¬ 
tions (black solid line) are in good agreement with the ex¬ 
perimental data (filled circles). Genuine non-Gaussianity is 
manifested above Sg = 8/3®^® (blue dashed line), (b) J for 
/|0)(0| + (1 — /)|2)(2| measured in experiment using the par¬ 
allelogram test with the fixed set of four phase-space points, ( 
-0.110, -0.110 ), ( 0.121, 0.100 ), ( 0.100, 0.121 ), and ( 0.331, 
0.331 ). Red circles (i), (ii), and (iii) are results for states 
|2)(2|, (|0)(0| + |2}(2|)/2, and |0)(0|, respectively, with their 
experimentally constructed Wigner functions shown below. 


state I I, n) to I t, + 1)) and a 7r-pulse for internal state 
transition | fi'n + 1) to | 4,,n -b 1). (h) We then dis¬ 
place the initial state |n = 0) or |n = 2). We perform 
a displacement operation T}{a) using the Hamiltonian 

Hd = -bh.c., realized by the same Raman laser 

beams with the beat frequency near wx and phase (j>, 
where the displacement is given by a = _ 

[ 2 ^ and ( 2 ^ for further description of experimental set¬ 
ting.) (iii) After each displacement D{a), we obtain the 
population Qn{o) = {n\D{a)pD'^ {a)\n) of Fock state |n) 
by observing | t)-state probability under the evolution 
with i?ajc and fitting the signal to [ 2 ^ 

(*) = ^ X! (1 “ Abe~^’’* cos (H„,„+it)) . (7) 

n 

^n,n+i ~ Irjfl is the Rabi oscillation rate be¬ 

tween I l,n) and | f; ^ + 1) while Ab and Xb are 
parameters incorporating experimental imperfections. 
We obtain the value of Wigner function ^W{x,y) by 
We observed the evolution with 
Hajc in steps of 1 ps for a total duration of 150 ys, about 
20 times the 7r-pulse duration of Hop, with 100 measure¬ 
ments at each step [ 2 ^. We repeated experiments 10 
times to obtain the parity under each displacement. 


We first perform an experiment on the state |n = 2) 
by changing the shape of parallelogram in phase space 
(Fig. 1 (a)), essentially the same as applying squeezing 
operation in Eq. The experimental data in Fig. 3 
(a) confirms the increasing capability of our test with 
squeezing to manifest genuine non-Gaussianity. 

In a second experiment, we obtained 77 = 2.20 ± 0.07 
for the ground state |n = 0) and J = 2.56 ± 0.1 for 
the state \n = 2), with the same set of four phase- 
space points [Fig. 3 (b)] [ 2 ^. A mixture of these states 
/|0)(0| -b (1 — /)|2)(2| can be addressed by including 
all data together for both states with an appropriate 
weighting [^. We observed 77 > Bq = 8/3®/® even 
for 0.5 < / < 0.65, where no negativity appears in the 
Wigner function (See for the experimental construc¬ 
tion of the Wigner function.) Due to the fluctuations 
of parameters in experimental system, we have standard 
deviation of parity about 5 % of the mean values, which 
could be further reduced by repeating experiments more. 

In summary, we developed a phase-space approach to 
manifest nonclassicality of a single-mode quantum state. 
Our approach has a close connection to the nonlocality 
test ruling out local hidden-variable models, particularly 
the BW-CHSH inequality, using that two quadratures 
of a coherent state behave like local variables in phase 
space. Our test was further extended to detection of 
genuine quantum non-Gaussianity, which is known to be 
an essential resource for many applications. Our methods 
successfully detect nonclassicality and non-Gaussinity for 
a broad class of Gaussian and non-Gaussian states even 
when the Wigner functions are positive definite, which 
was further confirmed by experimental data. Our phase 
space approach for single-mode fields and the BW ap¬ 
proach for multi-mode fields can therefore be crucial tools 
together to shed light on nonclassicality of GV systems. 


While we demonstrated our criterion using a trapped 
ion system, we anticipate that other CV systems can also 
be experimentally addressed, e.g. quantum optical sys¬ 
tems [3lLl32|. For example. Harder et al. recently showed 
an experimental method of obtaining photon-number dis¬ 
tribution robust against detection inefficiency by employ¬ 
ing only a thermal state as a probe and an on-off detector 
[32l[ . Combined with several displacement operations, it 
can be directly used to test our proposed criteria. 
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SUPPLEMENTAL MATERIAL 


A. TESTING GAUSSIAN STATES 

(i) Representing an arbitrary Gaussian state 

A Gaussian state a is fully characterized by its first-order 
moments (averages) ( 5 ) and (p) and second-order mo¬ 
ments (covariance matrix elements) 

+ R]Rk) — {Rj){Rk), j, k = 1,2 (8) 

where R = {q,p) with q = i(a’^ -I- a) and p = ^(d — 
d^). An arbitrary single-mode Gaussian state a = 
D{a)S{r, (j))ath{n)S^{r, cj))& {a) introduced in the main 
Letter has the averages (g) = Re [a], (p) = Im[Q!], and 
the covariance matrix elements 

Lii = i)(cosh 2r — sinh 2r cos 2(j)), 

^22 = ^ (d + ^) (cosh 2r + sinh 2 r cos 2(j)), 

ri 2 = r 2 i = —-(u-I-i) sinh2r sin2(/). (9) 


the displacement only shifts the center of the Wigner dis¬ 
tribution, while keeping its entire profile. We can thus 
ignore the first-order moments of a given state to find 
its maximum J'. Secondly, a phase-shift operation does 
not change the maximum either. If the maximum for 
p occurs at {xQ,yQ,xi,yi\ of four vertices of rectangle 
oriented at angle 9, the same maximum for the rotated 
state ap^-iipa^a gccurs at the vertices of the identi¬ 
cal rectangle oriented now at angle 9 p, because the 
phase-rotation also keeps the profile of Wigner function. 

(iii) M 2 Lxiniuni J for a Gaussian state 
Therefore, we set a = 0 to find a maximum JT without 
loss of generality. Introducing rescaled coordinates 

{x', y') = 2 p{x\/^ 2 , y (13) 

we recast Eq. (HH to 

Wa{x', y') = — exp[-a;'^ - -h kx'y'], (14) 

TT 

with 

J 2 ri 2 2 tanh 2 rsin 2(0 — (()) 

33 33 j 

riir 22 y 1 — tanh^ 2 r cos^ 2(0 — (()) 


Its Wigner function takes a Gaussian form as 

Wa{q,p) = — exp[-i(g - (g))'^r"yQ - (Q))], (10) 

TV Z 


where Q = {q,p)^, {Q) = (Re[a],Im[a])^, and p = 

If we rewrite Eq. m using the rotated quadratures 
Q = {x,y)'^ = 'R-{9){q,pY' (Fig. 1 (a) in the main Let¬ 
ter), we have 


Wa[x,y) = —eyip[-]-{Q - {Q)Y'T ^{Q-{Q))\, (U) 

TT Z 


The purity p = only affects the overall factor, with 
k independent of it. We thus set h = 0 for a maximum 
J'. For a given squeezing r, the parameter k is bounded 
by |A:| < 2 tanh 2 r, where the bounds are saturated at 
9 — (j) = ±-|, at which the rescaling in Eq. (13) reads 

ix',y') = p{x,y) with k = (the same k below 

Eq. (5) in the main Letter). 

We may be restricted only to /c > 0, as the maxi¬ 
mum for fc < 0 is equivalentely identified by y' ^ —y' in 
Eq. (14). Furthermore, as we show below, an optimal 
monotonically increases with k. This means that for a 
given squeezed state an optimal value occurs at 


where F = TZ{9)nZ{—9) represents the covariance matrix 
in a rotated frame, with 

R(9) = ('“* (12) 

^ ^ y— Sin 0 cos OJ ^ ^ 

and (Q) = (Re[ae“*®],Im[ae“*®]). Our task is then to 
optimize in Eq. (2) of the main Letter by varying the 
angle 9 in choosing a rotated rectangle and the coordi¬ 
nates {xo, poj aoi, 2 / 1 } of its four vertices. 

(ii) Invariance of maximal J under displace¬ 
ment and phase-shift operations 
We first note that the maximal value J for a given 
state p (both Gaussian and non-Gaussian) is the same 
as that of its displaced state D{a)p&{a): If the max¬ 
imum for p occurs at points {xo,yQ,xi,yi}, the same 
maximum then occurs for the displaced state at shifted 
points {x',y') = {x + Ke[ae~''^],y + because 


(16) 

i.e., the rectangle is rotated 45° from the squeezed axis. 
Intuitively, this choice of rectangle makes a sense because 
our test of nonclassicality relies on the non-factorizability 
of Wigner function. (See the paragraphs below Eq. (1) in 
the main Letter.) Eor instance, if we take 0 — (/) = 0, the 
Wigner function in Eq. (14) is factorized to Wfj{x', y') = 
W{x')W{y'), leading to no violation of Eq. (3) in the 
main Letter. 

—Condition for an optimal J 

Now the problem is reduced to the optimization of 

—g-Xg-vl+kxoyo _j_ ^-x^-yl+kxoyi 
pg-^l-vl+kxxyo _ ^-xl-yl+kxxyi 

where Q < k <2. The minimum Rk must be trivially — 1 
for all Gaussian states: Tk consists of four exponential 
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terms, each taking a value in [0,1] with — kxy = 

{x — + (1 ~ ^ 0. The minimum can occur at, 

e.g., ccq = yo “t oo and xi = yi = 0, which makes the 
first three terms zero and the last —1 in Tk- 

For a maximum, we obtain the partial derivatives 

^ = ^-^l-vl+kX0V0^_2x^ + fcyg) 
oxo 

+ e-^o-^i+'=“«^i(-2a;o + fcyi), 

BT, 22 

-^ = e-^o-yo+k-oyo^kxo-2yo) 
oyo 

+ e-^^^-Vo+^^^y°{kxi-2yo), 

^ + kyo) 

uXi 

^ = e-^o-yl+k^oyi _ 2y^) 

oyi 

- {kxi - 2yi). (18) 

Using the conditions = 0 at 

^ ® ^ oxq ayo oxi dyi 

optimal points, we have 


dH 

dxQ 

dTk 

dyo 

m 

dxi 

dTk 

dyi 


{kxo - 2yo){-2xi + kyo){kxi - 2yi) 
(-2x0 + kyo)(-2xi + kyo){kxi - 2yi) 
(-2x0 + kyo)(kxi - 2yo)(kxi - 2yi) 
(-2x0 + kyo)(kxi - 2yo)(-2xi + fcyi) 


=2fc(4 - e)(xo - xi)(yo - yi)(xoXi - 

=0. (19) 


y 



where f(c,t,k) = Obviously, if one finds a 

global maximum for in Eq. (22) with a set of values 
{y',c, t, fc}, the same maximum for in Eq. (21) must 
occur at {y = -^^=^^,c,t,k}. Now, taking a partial 
derivative of T over k, 


dk 


2y' ^f 02g-tV(c,t.fc)y'^ _^g-/(c.i.fc)y'^'\ ^ 

(u +1 — kty \ j 

(23) 


we see that J-k is a monotonic function of fc € [ 0 , 2 ] for 
any given values of {y',c, t}. This means that we only 
need to consider either fc = 0 or fc = 2 at each point 
of parameter space to find a global maximum. That 
is, max{J^fc=o(y^ c, t), .Ffe= 2 (y^ c, t)} becomes the global 
maximum of T by covering the whole ranges of {y', c, t\. 
As fc = 0 is the case of coherent state yielding only a 
classical bound 2 , we conclude that the case of fc = 2 , i.e. 
an infinitely squeezed state, gives the maximum bound 
for all Gaussian states. 

—Maximum Gaussian bound 
Coming back to Eq. (18) with k = 2, the partial deriva¬ 
tives are rather simplified to 


As Xg = Xi and yo = yi give only a classical value, Tk = 
2e~^o-yo+k^oyo < 2, the optimal points we seek must 
satisfy xoxi = yoyi, which is equivalent to 


^ = EL =f 
yo xi 


( 20 ) 


thus Xo = tyo and yi = txi. Furthermore, we may set 
xi = cyo, which enables us to cover all pairs of {xi,yo} 
by scanning c in the entire range of real values. 

Then, setting (xo,yo,Xi,yi) = (ty,y,cy,tcy), in 
Eq. (17) reads 


J-fe(y, c, t) = 

By defining y' = y\/y + \ — kt, the above becomes 
Tk(y', c, t) = e-y'^ - -f -b 


^ = -2e-("'>-^“)'(xo - yo) - 2e-^-°-y^^\xo - yi), 
axo 

= 2e~^^°~y°^^ (xo -yo)+ (xi - yo), 

oyo 

^ = -2e-("^-^'>)'(xi - yo) + 2e-("i-^^)'(xi - yi), 

UXi 

^ = 2e-("“-^^)'(xo - yi) - 2e-("i-^^)'(a;i - yi). 

oyi 

(24) 


A trivial solution for = 0 is 

Xo = yo = Xi = yi which only yields T 2 = 2. To find the 
other solution, we make two observations. (1) The above 
equations impose an ordering for optimal parameters. If 
we have xi > yi or xg < yo, substituting the condition 
to the equations yield xi > yo > xo > yi or xi < yo < 
Xo < yi, respectively. From now on, we assume Xi > yi 
, without loss of generality. (2) Every equation is in the 
(22]^ame form; ae““ —/3e~^ = 0. From = 0, we 
























obtain x^ — yi = xi — t/o in view of the ordering of optimal 
parameters and Fig. |4l As can be seen from Fig. |4l the 
number of solutions for xe~^ = 5 is at most two. In 
addition, none of Xq — yi and Xi — y^ can be identical to 
xi — yi because of the ordering. The only possible case 
is then x^ — yi = x\ — y^y = a and xi — j/i = /3 where 
a < P with Q;e““ = /3e“^ . Using a similiar argument 
to ^ = 0, we also have Xo — yo = —a. Combining 

the results, we have (3 = 3a; xi — yi = (xq — yi) + (xi — 
yo) — (xq — yo)- It enables us to rewrite the function as 
= 3e““ — , which gives the maximum at 

a = Combined with the condition XqXi = yoyi, 

it also gives the optimal points xg = —yo = —^o; and 
xi = -yi = §a. 

—Optimal for finitely squeezed states 

Let us now find an optimal Tk in Eq. (21) for each k £ 
[0, 2]. Its partial derivatives are given by 


+ 2t(c2 - fcc+ 


dTk 

dc 




^ = -2y{if -kt + 
dy 


+ (c 2 -fcc+l)e-("'-'="+^)^' + 
- if -kt + 


(25) 


In the case of fc = 2 considered before, the optimal points 
were given as xq = —yo = —^ 0 ( and xi = —yi = |a. 
This indicates t = —1, recalling the parametrization 
{xo,yo,Xi,yi) = {ty,y,cy,tcy). Indeed, we find that a 
set of optimal points satisfy = 0, with 

t = — 1 and 


g-{(l+fe)c^+feo-l}y^ 


g-(l+c)(l-c+fc)y= 


2c — k 
c(fc + 2)’ 
kc — 2 
k + 2' 


which gives 


p _ 0 

{(l + fc)c2 + fcc-l}l0g^^ 

= (lTc)(l-cT.)log^ 


(26) 


(27) 


For a given k £ [0, 2], we can numerically obtain c from 
Eq. (27), and then y from Eq. (26). (The parameters Zi 


introduced below Eq. (5) in the main Letter then reads 
zq = y and Zi = cy.) The resulting maximum Tk (blue 
solid curve) is shown in Fig. 2 (a) of the main Letter, 
which is also in agreement with the result obtained inde¬ 
pendently from a full numerical optimaization (red dots) 
by scanning all real ranges of {t, c, y} in Eq. (21). 

As an analytical example, Eq. (27) gives c = 3 for 
fc = 2, which recovers the optimal points xg = —yo = 
— ^a and xi = —yi = |a obtained independently in the 
previous subsection. For a mixed state with y < 1, we 
only need to multiply a factor y = as seen from 

Eq. (14), to a result for a pure state with the same r. 


B. RIGHT TRIANGLE TEST 

We here introduce another nonclassicality test taking 
three phase-space points into consideration, which has no 
analogue in nonlocality test, as 

J'[p\ = |{IF;(xi,yo) + <(xo, 2 /i)-<(xi, 2 /i)}, (28) 

where the three points may form a triangle rotated by 
an angle 0 from g-axis in phase space. Our test is linear 
with respect to a convex mixture of quantum states, i.e. 

displacement does not 
change the optimal value of like the four-points test in 
the previous section. Investigating a vacuum state thus 
suffices to find its classical bound as 

J7'[|0)(0|] = -f (29) 

Regardless of the sign of we find the clas¬ 

sical bound as J'[pc\\ < 1; 

(i) if e- 2 “o - e- 2 ^? < 0 , we have J^'[| 0 )( 0 |] < 
^-2xl^-2yl ^ 

(ii) if > 0 , we have 

:^'[| 0 )( 0 |] < e-^^"e-^y- + e- 2 "o - 

= _e-2^?(l-e-2!^o)-he-2^o < 1. (30) 

In contrast, there does not exist a useful nonlocality 
test constructed out of a linear sum of three product val¬ 
ues only, as B' = aibg -|- agbi — aibi. One can readily 
find a set of local values that achieve a maximum value 
3 (oq = l,ai = —1 and bg = —l,&i = 1 ) and a mini¬ 
mum value -3 (ag = —l,ai = 1 and bg = —l,bi = 1 ), 
respectively. Therefore, our proposed test here taking 
only three points is a unique nonclassicality test having 
no analogue in nonlocality test. 

For a given Gaussian state with (a,r, 0, h), the maxi¬ 
mal value for this triangle test occurs at, 

Xj = a^: + (-lyz'K 

Vj =Oiy- (-I)^z'k, 


( 31 ) 
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where j € {0,1}, 0 = <!> + ^, k = and Zj can 

be found analytically, similar to four-points test. 

—Optimal J'' 

Here, after the rescaling of coordinates ^{x,y) —>■ {x,y), 
we need to optimize 

_ ^-^l-vl+kxoVi _|_ ^-xl-vl+kx-iya _ ^-xl-yl+kxiyi 

The conditions ^ ^ ^ ^ = 0 im¬ 

mediately give xi = Ij/o and yi = |xo. Then, us¬ 
ing a parametrization {xQ,yo,xi,yi) = {ty,y,‘ly,\ty), 
we consider 'Hk{y,t) = -p — 

e~(* ^ which leads to = 0 as 

, k'^, 2t , , 2 ,, 

The combination of the above two equations gives 
(kt - 2r log[(l - ^)^] = ( 2 t - log[(l - 

(34) 

Noting that 'Hkiy, —Kl) > 'kLk{y, |t|), an optimal point 
must occur for t < 0. We have numerically confirmed 
that < = — 1 is the only negative solution to Eq. (34) for 
any k G [0,2]. Then, we obtain y^ = , log[—and 

1 2 —fc 

the maximum'Hfc = — ^(2 —fc) 2 +v ( 2 -|-/c) monotonically 

increasing to 2 with fc. The optimizing parameters in Eq. 
(31) read z'q = y and z[ = ^y. 

Similar to the rectangle test using four phase-space 
points, the triangle test detects every pure Gaussian non- 
classical state. Moreover, Eig. 2 (b) of the main Letter 
clearly demonstrates that the triangle test outperforms 
the rectangle test to detect nonclassical Gaussian states. 

Using the fact that no Gaussian state overcomes the 
bound Bq = 2 , we can also construct a genuine quantum 
non-Gaussianity test as 

yv' = j'[UGpuU-^, (35) 

where a positive value of the witness, W' > 0 for a cer¬ 
tain Ug, manifests that the state is genuinely quantum 
non-Gaussian. In order to beat this bound, Wjy^^jy-t {q,p) 
must necessarily take a negative value. While this test 
may be interesting in itself, it would not be practically 
significant compared with the rectangle or the parallel¬ 
ogram test in the main Letter as the latters can detect 
quantum non-Gaussianity of a positive Wigner function. 


C. VIOLATION OF LOWER BOUND -1 FOR A 
NONCLASSICALITY TEST 

The lower bound —1 of our nonclassicality test, Eq. (3) 
in the main Letter, can be more useful than the upper 



FIG. 5: A noisy single photon state, /|0)(0| -f (1 — /)|1)(1| 
under rectangle (blue solid) and parallelogram (red dashed) 
tests. The violation of lower bounds demonstrates nonclassi¬ 
cality and non-Gaussianity better than that of upper bounds 
(inset). That is, the lower bounds detect the state in the 
ranges 0 < / < 0.25 and 0 < / < 0.33 under rectangle and 
the parallelogram tests, respectively, while the upper bounds 
do so in 0 < / < 0.08 and 0 < / < 0.01. 


bound 2 if there exists a substantial negativity in the 
Wigner function of the original state, even when it is 
mixed with some noise. For example, we plot the results 
for a noisy single photon state, /| 0 )( 0 | + (1 — /)| 1 )( 1 |, in 
Fig. [S] It shows that the violation of the lower bound 
occurs for a broader range of / than the violation of the 
upper bound. 


D. MAXIMAL VALUE OF J[p] IN EQ. (2) BY 
QUANTUM NON-GAUSSIAN STATES 


In order to find out a maximum possible value of by 
a quantum non-Gaussian state in Eq. (2), we here solve 
an eigenvalue problem = X\ip), where the Hermitian 
operator H is defined as 

1 1 

n = ^ ^(-l)“'U(g, + ipt){-lfDHqa + ipb). (36) 

a =0 6=0 


We note that the four terms in T-L contributing to a linear 
sum have a different structure from those in the Bell op¬ 
erator to test nonlocality, B = dobg + dobi + Si&o — ai&i. 
In the latter case, using d| = = / {i,j = 0,1) 

and the commutativity of spatially separated observables 
[di,bj] = 0 , we are led to B'^ = AI + [do,ai] G) [5i,5o], 
from which the maximum (H^)max = 8 known as the 
Girelson bound 23[ can be derived. In contrast, each 
term of H in our test is not factorized to a product of 
two commuting observables (unlike the form dibj in B), 
which instead involves local non-commuting observables 
q,p and n. This indeed makes it possible to achieve a 
maximum eigenvalue of T-L above 2 -\/ 2 , as we show below. 

We may set {qQ,po,qi,pi) to be {0,0, d,d) without loss 
of generality. Note that we can always transform arbi¬ 
trary four variables {qo,po,qi,pi) to (0, 0,(i, d) by Gaus- 
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(a) Anti-Jaynes-Cummings Coupling (b) Displacement Operation 



FIG. 7: Experimental configuration with external driving 
fields to produce (a) the anti-Jaynes-Cummings interaction 
between internal and motional degrees of freedom and (b) 
the displacement operation on the motional mode. 


which is shown in Fig. [6] (b). 

As an illustration, we show the case of = 1, i.e. a 
superposition of 3 x 3 coherent states in Fig. |6] (c) and 
(d), which achieves a value of J[p\ ~ 3.32. 


FIG. 6: (a) Maximum eigenvalue Ajv of H in Eq. (36) that 
is numerically obtained by truncating n and m up to A in 
the recurrence relations of Eq. (38). (b) Expectation value 
/Tiv = {H) of a truncated superposition of (2A +1) x (2A +1) 
coherent states in Eq. (I39II compared with Ajv. (c) Wigner 
function and (d) its contour plot for a superposition of 3 x 3 
coherent states with a lattice size 2d — in Eq. (I39II . 

sian unitary operations, which do not affect the maxi¬ 
mum value of J achieved by a non-Gaussian state. We 
seek a solution in a form 

OO OO 

i^)= E E C'„,m|2d(n-b fm)), (37) 

n=—OO m——oo 

where |2d(n -I- im)) is a coherent state (n and m: in¬ 
tegers), that is, a superposition of coherent states in a 
2-dimensional lattice of size 2d. Using Eqs. dMl) and dSil) 
and the identity i)(a)(—l)”Dl(a)| 7 ) = e““^*+“*'^| 2 a — 
7 ), we obtain a recurrence relation 

^ 1 —Aid^m^ I 4:id^n^ 

—^ — n, — m “r c O — — “h o 0_yi^_77i,-(_i 

-Aid^{m-n) ^ /qo\ 

In Fig. | 6 ](a), we show the maximum eigenvalues numeri¬ 
cally obtained by a truncation of numbers up to N in the 
recurrence relations, which turns out to approach 4 with 
N increasing for a choice oi = ^ + j {R: integer). 
To check whether these values can actually be achieved, 
we compare Xn and where 

N N 

1 '*/’^) = E E Cn,m\2d{n + im)) , (39) 

n— — N m— — N 


E. EXPERIMENTAL REALIZATION 

The control and manipulation of phonon mode can 
be achieved by the anti-Jaynes-Cumming interaction 
that couples the phonon mode d to the internal states 
represented by Pauli operators <7^, which is given by 
i7ajc = ^d^(7+-|-h.c.. It can be realized by two counter- 
propagating Raman laser beams that have the beat fre¬ 
quency near whf + wx as shown in Fig. [3 (a). On the 
other hand, for a displacement operation on the phonon 
mode, the Hamiltonian Hd = cd -bh.c. can be re¬ 

alized in the same experimental configuration, with the 
beat frequency near wx and the phase as shown in Fig. 

[21(b) [ 23 . 

For the measurement of J and Wigner function, first 
we perform a displacement operation on an initial state, 
then apply the anti-Jaynes-Cumming interaction to ob¬ 
tain the phonon number distributions, which provide the 
value of the parity, Figure [5] shows typical 

time evolutions of Pp), i.e., population in the state | f) 
under the operation of anti-Jaynes-Cummings interaction 
Hajc and the phonon-number distributions Qn that are 
obtained by fitting the time evolution signal to Eq. (7) 
in the main Letter. 

We also experimentally measure the density matrices 
of the state |0)(0| and |2)(2| and construct their Wigner 
functions using the density matrices. To measure the 
density matrix, we basically follow the method in Ref. 
281 and apply the iterative maximum-likelihood method 
33j to avoid unphysical, non-positive, matrix by the di¬ 
rect reconstruction. In experiment, we displace the state 
p to eight different angles with the fixed amplitude of 
a = 0 . 8 . 
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FIG. 8: (a)(c) Typical time evolutions with the application 

of Hajc and (b) (d) phonon distributions after fitting the sig¬ 
nal to Eq. (7) in the main Letter. (a)(b) are the results 
after displacement operation on a vacuum state |n = 0) to 
(0.331,0.331) in phase space, while (c)(d) are the ones after 
the same displacement operation on the state |n = 2). 



FIG. 9: Contour plot of Wigner function for a Fock state 
|2), with (a) four points chosen for a nonclassicality test, and 
(b) four points (green) of rectangle transformed to those (red) 
of parallelogram. Red points were used in our experiment of 
non-Gaussianity test. 


—Choice of phase-space points for test 


For our test of non-Gaussian states, there would not be 
a unique, systematic, way to choose optimal phase-space 
points, similar to the situation with nonlocality tests. 
However, we may give a rough description of choosing 
appropriate points as follows. For a rectangle test, as 
three terms contribute positively and one term negatively 
in J^ we may configure the rectangle such that one vertex 
corresponds to a negative value of Wigner function. An 
example of Fock state |2) is illustrated in Fig. IH] (a), 
together with the choice of four phase-space points. The 
left-upper point is placed at a negative Wigner function. 

On the other hand, when a pure non-Gaussian state is 
mixed with noise, e.g. /|0)(0| + (1 — /)|2)(2| in our ex¬ 
periment, the parallelogram test provides a much more 
efficient method. We may first choose a rectangle with a 
modest size (green dots in Fig. |9] (b) ) and then shift the 
points (squeeze the rectangle) according to the rule given 
below Eq. (6) in the main Letter, resulting in a parallel¬ 
ogram (red dots in Figl9](b)). Note that one red point 
now reaches a negative region of |2) by this transforma¬ 
tion, although it is degraded by the vacuum state |0) (not 
shown) that has a positive-definite Wigner function. Of 
course, we need optimization over the initial rectangles 
followed by squeezing, but it can be done at least nu¬ 
merically with some prior information on the produced 
state. This is how we chose the phase-space points for 
/|0)(0| -I- (1 — /)|2)(2|. Even when selected points are not 
strictly optimal, we can experimentally demonstrate the 
genuine non-Gaussianity to a large extent, as the inequal¬ 
ity (6) of the main Letter is satisfied in a rather broad 
parameter region. In fact, this was the case with our ex¬ 
periment, taking four (non-optimal) phase-space points 
well-separated for an easy control of displacement. 













